%Bedoeling van dit script is voor een bepaald net (met een vaststaand aantal en soort sectoren per bus) de kosten tgv PQ
%problemen te berekenen PER BUS.
%Deze gegevens worden opgeslagen in een array Kosten_Per_Belastingsbus. In
%de eerste kolom hiervan staat het nummer van de belastingsbus, in de
%tweede kolom staat de PQ kost op die bus. Na het runnen van
%kosten_mitigatiemiddelen_net_deel2.m staat in de derde kolom de kost voor
%mitigatiemiddelen en in de vierde kolom de ratio van de baten naar de
%originele kosten.

load belastingsbussen  %gegevens over de verschillende belastingsbussen laden
load problemen
sectoren_namen; %namen sectoren laden
[NUMERIC1,TXT,IndelingBussen] = xlsread('Verdeling_belastingsbussen.xls');
[r,k]=size(IndelingBussen);
kost_deze_bus=0;
totale_kost=0;
Kosten_Per_Belastingsbus = zeros(aantal_belastingsbussen,4);
for i=1:1:aantal_belastingsbussen
   nummer_belastingsbus = belastingsbussen(i,1);
   Kosten_Per_Belastingsbus(i,1) = nummer_belastingsbus;
   rijen=find(cell2mat(IndelingBussen(2:size(IndelingBussen,1),1))==nummer_belastingsbus)
   nummer_belastingsbusstr = int2str(nummer_belastingsbus);
   aantal_sectoren_op_deze_bus=size(rijen,1)
   for l=1:1:aantal_sectoren_op_deze_bus
       rij=rijen(l,1);
       naam_sector=char(IndelingBussen(rij+1,3));%Plus 1 want in de eerste rij van het excel sheet staan titels
       test_string=strcat('kosten_sector_',naam_sector);
       aandeel_deze_sector=cell2mat(IndelingBussen(rij+1,4))/100;
       temp=char(test_string);
       test=eval(temp);
       kost_deze_sector=test(i,2)*aandeel_deze_sector;
       kost_deze_bus=kost_deze_bus + kost_deze_sector;
   end
   Kosten_Per_Belastingsbus(i,2) = kost_deze_bus;
   totale_kost=totale_kost + kost_deze_bus;
   kost_deze_bus=0;
end
